Magnetic Raman Scattering of Insulating Cuprates 
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We study the B lg and Ai g Raman profiles of M2CUO4 ( with M= La, Pr, Nd, Sm, Gd), 
Bi2Sr2Cao.5Yo.5Cu208+;,, YBa2Cu306.2 and PrBa2Cu2.7Alo.3O7 insulating cuprates within the 
Loudon-Fleury theory, in the framework of an extended Hubbard model for moderate on-site 
Coulomb interaction U. We calculate the non-resonant contribution to these Raman profiles by 
using exact diagonalization techniques and analyze two types of contributing mechanisms to the 
line shapes: 4-spin cyclic exchange and spin-phonon interactions. Although these interactions con- 
tribute to different parts of the spectra, together, they account for the enhanced linewidth and 
asymmetry of the B\ g mode, as well as the non-negligible intensity of the A\ g Raman line observed 
in these materials. 

Pacs Numbers: 78.30.-j, 71.10.-w 



I. INTRODUCTION 

Inelastic Raman and neutron scattering experiments 
on insulating cuprates provide information about the spin 
dynamic of the CuC>2 planes of high-T c precursors [jE|. 
The observed short and long wave-length low-energy ex- 
citations have been described using the two-dimensional 
spin-i antiferromagnetic Heisenberg model (AFH), and 
the standard theory of Raman scattering based on the 
Loudon-Fleury Q (LF) coupling between the light and 
the spin system. Within this theoretical scheme, the 
value of the Cu-Cu exchange constant J has been es- 
timated from the first moment M\ of the B\ g Raman 
line and, in virtual agreement, from neutron scattering 
measurements of the spin velocity. For insulating ma- 
terials having the Cu atom on different chemical envi- 
ronments, the estimated value of J results in the range 
- (0.10-0.14) eV. 

Within the AFH framework, numerical and analytical 
calculations have been able to describe, the temperature 
dependence of the spin-spin correlation length obtained 
from neutron scattering data (|,||], the temperature de- 
pendence of the spin lattice relaxation rate 1/T± mea- 
sured by Nuclear Magnetic Resonance , and a promi- 
nent structure, ascribed to 2-magnon scattering, of the 
Big Raman line j7[-pd[|. Although, all these theoretical 
results suggest the AFH on the square lattice as a very 
good starting point for describing the spin dynamic of the 
insulating cuprates, some features of the Raman lines re- 
main to be explained within the LF theory. In particular, 
it is not yet understood the AFH-LF failure to describe 
the Bi g spectral shape and its enhanced width. And of 
course, the AFH-LF can not reproduce a non- vanishing 
A\ g response. 

The shape of the Bi g line presents some characteristic 
features common to several insulating cuprates. In fact, 



at room temperature, a very similar line shape has been 
observed for all members of the M2CUO4 series (M= La, 
Pr, Nd, Sm, Gd ), Bi 2 Sr 2 Ca .5Yo.5Cu 2 8+y (BSCYCO), 
YBa 2 Cu 3 6 . 2 and PrBa2Cu3.7Alo.3O7 (PRBACUALO) 
|l2| p^| . In all cases, the line extends up to ~ 8 J, which 
in the Ising limit, can be ascribed to multimagnon excita- 
tions involving 16 spins. Unfortunately, calculations on 
the AFH show a negligible contribution of multimagnon 
processes not only for the Raman line but also for the 
spin structure factor at the antiferromagnetic wave vec- 
tor. Another feature is a very intensive peak located at 
lo ~ 3 J, identified as 2 — magnon scattering ]^-|l0|]. Now, 
it is believed that the asymmetry of the line originates 
on a second hidden peak on the high energy side of the 
spectrum at uj ~ 4J whose intensity is ~ 25% smaller 
than the intensity of the 2-magnon peak [Q. In the 
Ising limit, the energy of these magnetic excitations cor- 
responds to 4-magnon scattering. For the AFH model, 
they give a very small contribution to the intensity of 
the spectrum. On other geometries interesting features 
also appear. The Ai g line has a maximum at a slightly 
high frequency with respect to the 2 — magnon peak [ fL7| , 
while for the Ai g and B^g symmetries the center of the 
spectrum is located at lu ~ 5J and u) ~ 4J respectively 
ifLsf . Finally, it has been found a slight temperature de- 
pendence of the scattering intensities of the Bi g line in 
the insulating cuprates. In fact, the rise in temperature 
from 30 to 273K causes the 2 — magnon peak intensity 
to decrease by only ~ 10%, in contrast to ~ 50% for the 
5 = 1 system La2Ni04 jl2f - Consistently, previous calcu- 
lations show a weak temperature dependence of the B\ g 
spectrum for T < J |lj]] for the spin | AFH model, which 
in turns suggest the importance of short-range spin-spin 
correlations. 

In the Fleury-Loudon theory, the state emerging from 
the application of the current operator on the ground 
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state is considered as an eigenstate. In this intermediate 
eigenstate, the electron and the hole, produced by the 
photon scattering, are very close in real space, creating a 
charge transfer exiton, and the individual propagation of 
each particle is neglected. This approximation could be 
appropriated for the case of the cuprates due to the pres- 
ence of an electron-electron interaction (U p d) between 
nearest neighbor sites (oxygen and copper sites). This 
interaction, not included in the Hubbard Hamiltonian, 
is important to describe the charge excitations of the 
cuprates and induces an attractive interaction between 
the hole and the electron | p0| . Due to this interaction, 
the hole and the electron are close in space and there is 
not free individual propagation of them [|21|. Further- 
more, the disorder present in the cuprates [|22[ tends to 
localize the electron-hole pairs. Chubukov and Frenkel 
tried to include the resonant effects in the frame of the 
one band Hubbard model [^3|. Schonfeld et al. |24j] cal- 
culated the resonant contributions but they could not 
properly describe the B\ g profile, and the observed A\ g 
line remains to be explained. In both cases, they did not 
include the effect of the U p d interaction to describe the 
charge excitations and considered that the electron and 
the hole propagate individually. For the reasons given 
above, we think that the Fleury-Loudon theory is an ap- 
propriated framework to describe the Raman spectra of 
the cuprates. 

Even when the resonant effects may be needed to ex- 
plain the dependence of the Raman spectra on the incom- 
ing laser frequency, the main features of the line studied 
in our work persist over the wide range of frequencies 
used in the experiments (l^JT^,|I^,|l^,^5|,^6| and are es- 
sential the same than out of resonance. This point is 
shown in, e.g., the recent work in PrBa 2 Cu2.7Al .3 |ll| 
where the Raman spectra out of resonance was measured 
and the main features of the B\ g line shape remain un- 
altered. Besides, the LF theory has been the standard 
framework to describe the Raman spectra of other an- 
tiferromagnetic S = 1/2 compounds [^,^8|. So, it is 
an interesting issue whether the persistent aspects of the 
Raman line can be well described within the LF theory 
with an appropriate model. 

Previous theoretical work using LF theory tackled the 
problem of the line's shape by using series expansions, 
variational Monte Carlo, interacting spin-waves and ex- 
act diagonalization techniques on small clusters |t] [ll]] . 

Different theoretical scenarios have been proposed for 
describing the anomalously enhanced width of the B\ g 
line, and although at a first sight, quantum fluctuations 
could be the main contribution to the width of the B\ g 
line, this last hypothesis has been rejected after mea- 
surements on the S = 1 system MPS3 for which a rel- 
ative width comparable to the insulating cuprates has 
been observed p9| . A different proposal has been moti- 
vated by the temperature dependence of the 2-magnon 
peak, namely, an increase of the Raman linewidth with 



increasing temperatures [B0|. This has been considered 
as a strong indication of a phonon mechanism for the 
line's broadening and its effect has been analyzed through 
a non- uniform renormalization of the exchange constant 
J pjj j. Another argument which support the phonon 
mechanism is that although the second cumulant Mi 
of the Bi g Raman spectrum is almost J- independent 
for the materials mentioned above, its exchange constant 
changes by ~ 20%. This indicates a linewidth dominated 
by non-magnetic contributions such as intrinsic disorder 
of the spin lattice, defects or phonons. Of course, this 
does not rule out other magnetic mechanism providing 
scattering at frequencies higher than the 2-magnon peak. 

So far, the phonon mechanism alone is insufficient to 
describe the line asymmetry due to important contribu- 
tions of 4-magnon scattering Eh}. This requires to go 
beyond the minimal AFH model. An effective descrip- 
tion of the insulating cuprates based on the single-band 
Hubbard model p2| , recently appeared, gives a theoreti- 
cal framework to go beyond the AFH model and provides 
additional multispin interactions which finally could con- 
tribute not only to the Raman spectra but also to the 
midinfrared phonon-assisted optical conductivity ]35} | . In 
this work, we extend previous calculations of the B\ g and 
A\ g lines based on the minimal AFH model by includ- 
ing multispin and spin-phonon interactions. We calculate 
the Raman spectra on finite-size systems using the now 
standard Lanczos method |34) for the series of materials 
mentioned above. Since inelastic light-scattering mea- 
surements probe short-wavelength spin fluctuations, we 
expect finite size effects to be small p5[ ]. We find that 
two mechanisms, multispin and spin-phonon interactions, 
contribute to the width and shape of the B\ g Raman line, 
and at the same time, allow an otherwise forbidden A\ g 
response. 

II. RAMAN SPECTRA AND EFFECTIVE SPIN 
MODEL 

The scattering of light from insulating antiferromag- 
nets at energies smaller than the charge-transfer gap 
~ (1.5 — 2)eV can be described phenomenologically, by 
introducing a Raman operator, based on the symmetries 
of the magnetic problem. For the one-dimensional irre- 
ducible representations of the square, these operators are 
given by 

^ B is = X/ ^ i '('^' i + e * - 
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°A 2g = y^e AW 5 i .(5i + e ( , x S- 1+ev ). 

i 

where /Lt = ±x,±j/, and e M „ = — e„ M = -e_ M „. Here e x 
and e y denote unit vectors in a;,y directions of a square 
lattice. The Raman Hamiltonian is proportional to one 
of these operators, being the prefactor matrix elements of 
the dipolar moment and therefore depend on microscopic 
details of the system. Here, we will assume all them equal 
to one and therefore we will be unable to compare the rel- 
ative intensities of spectra of different symmetries. Mi- 
croscopically, these Raman operators can be derived from 
the theory of Raman scattering in Mott-Hubbard systems 
p6| . In the strong-coupling limit of the Hubbard model, 
they appear quite naturally as the leading-order contri- 
butions in an expansion in A = t/(U — u>). Of course, for 
the 2D AFH model, the Oa x operator commutes with 
the Hamiltonian, so the A\ g response vanishes at the low- 
est order and terms of order A 3 are required to have a 
small but finite response. 

The Raman spectrum Rr(u>) at T = 0, is given by 

R r (u) = J2\(MOr\^n)\ 2 S(iu~E n + E ) 

n 

where T is Ai g , Bi g , Ai g or Bi g , 1 0o > is the ground- 
state, Eg is the ground-state energy, and \<f> n > is an 
excited state of the system with energy E n . The spec- 
trum R^{u>) can be obtained from a continued fraction 
expansion of the diagonal matrix clement of the resolvent 
operator l/(uj + E + iS — H) between the state Or\4>o >■ 
Here S is a small imaginary part added to move the poles 
away from the real axis J34|. 

Let us describe first the B\ g Raman spectrum of the 
2D antiferromagnetic Heisenberg model. It has mainly 
two contributions: (a) a very strong two-magnon struc- 
ture located at uj ~ 3J and (b) weaker J^-magnon pro- 
cesses at higher-energies. In a recent work, A. Sandvik et 
al. calculated numerically the Heisenberg B\ g spec- 
trum for lattices up to 6x6 sites. They found a change 
in the strong structure at u ~ 3J in the exact spec- 
trum going from 4x4, with a single two magnon domi- 
nant peak, to 6x6, where two equal sized peaks appear. 
Unfortunately this is the bigger system available nowa- 
days and it is not clear if these two peaks are not an 
spurious characteristic of this cluster since this kind of 
structure is not observed in bigger systems calculated 
with Monte Carlo-Max Entropy jtl). They concluded 
that the Heisenberg model can not describe the broad 
Bi g spectrum and of course it gives a zero A lg pro- 
file even for that size of clusters. The 4 — spin exci- 
tations are generated by flipping two pairs of spins of 
the Neel background in a plaquete or in a column with 
an energy cost of 4J and 5 J respectively. The intensity 
of these high-frequency processes is very small < 10% 
P,B5|, and cannot account for the asymmetry observed 



for this mode |3jj] on the high frequency side of the spec- 
trum. Since a distinct shoulder appears at w ~ 4J, it 
has been proposed that the 4 — magnon channel is en- 
hanced by 4-spin cyclic exchange interactions ]l2|] . In 
fact, finite cluster calculations (3| |4(| for the undoped 
compound have shown that although within the Cu02 
planes the leading term in the effective spin Hamiltonian 
is the antiferromagnetic Heisenberg exchange interaction, 
other terms are required to describe the low energy spin 
excitations. These higher-order terms appear naturally 
by performing a canonical transformation up to fourth 
order on the Hubbard model jll|. A single-band Hub- 
bard model can simultaneously describe the low-energy 
charge and spin responses of insulating S^CuG^C^ |p2| . 
Hence, our starting model to describe the spin dynamic 
of all the compounds mentioned before, is an extended 
Hubbard Hamiltonian, including next-nearest-neighbor 
(NNN) hopping t' . By means of a canonical transfor- 
mation to fourth order in t and second order in t', the 
effective spin Hamiltonian is written as follows, 




+K (SiSjXSfcS,) + (SiSiXSjSk) 

-(SiSfcXSj-so 

with J = At 2 /U - 64i 4 /[/ 3 , J' = U l2 /U + U 4 /U 3 , J" = 
4i 4 /t/ 3 , and K = 80i 4 /[/ 3 . 5 5' 5" runs over NN; 2nd NN; 
3th NN and (i, j, k, 1} means the sum over groups of four 
spins in a unit square. The effect of the cyclic exchange 
interaction on the antiferromagnetic ground state of the 
2D Heisenberg model has been studied showing that for 
small values of K/ J, the antiferromagnetic phase is sta- 
ble JlU . The staggered magnetization increases as K/J is 
increased, has a maximum at ~ 0.75, and then decreases 
showing a transition to a spin canted region. In partic- 
ular, for the value of K estimated for the cuprates, the 
system is still in the antiferromagnetic phase and its stag- 
gered magnetization is m st ~ 0.58. The reduction of the 
three-band model onto the single band Hubbard model 
for realistic values of the multiband parameters, indicates 
that the effective hopping t is bounded between 0.3eV 
and 0.5eV while U/t ~ 7 - 10 Q. Furthermore, the 
derivation of the one band Hubbard Hamiltonian given 
by Simon and Aligia ( see Ref. |Q ) for the parame- 
ters obtained from LDA calculations for La2Cu04, gives 
t - 0.45ey, U/t - 7.6, t'/t - 0.15. Deviations from these 
values are expected for different materials, in particular, 
for the values of the long-range hoppings. In Table I, 
we show the values of the parameters for the 7 systems 
studied in this work. 
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III. RESULTS 

Let us examine first the role of multispin interactions 
on the spectra. Figure 1 shows the calculated Raman 
spectrum at different symmetries for U = 9.M and t' — 
0.28i. The B\ g line shows the same features as for the 
Heisenberg model, being the main difference the gain of 
spectral weight at high-frequencies, Fig. 1(a). The 4 — 
magnon processes are now enhanced having a relative 
intensity to the 2 — magnon peak of ~ 20%. This is in 
agreement with Ref. where it has been found that 
the deviation from the symmetric part of the B\ g mode is 
mainly due to a peak located at to ~ 4 J whose intensity 
is about 20% of the 2-magnon peak. For the A 2g and 
B 2g modes, Fig. 1(c), the center of gravity of the spectra 
agree rather well with the experimental data Eq] . 
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FIG. 1. Big and A\ g Raman spectra for 20 and 16 clus- 
ters. Calculations were done for the effective spin model with 
no spin-phonon interaction. The microscopic parameters cor- 
responds to YBa2Cu306.2 ( see Table I ). 

However, calculations for the Ai g line, Fig. 1(b), shows 
that the maximum is shifted by ~ J respect to the two 
magnon peak, in disagreement with the experimental 
shift. Although the 4-spin cyclic exchange interaction 
provides a mechanism for some features of the B\ g mode 
and yields responses in the correct frequency range for the 
A 2g and B 2g lines, it does not give correctly the frequency 
range for the Ai g spectrum neither the full linewidth of 
the Big line. As was mentioned before, this can be due 
to sample inhomogeneities, additional strong interactions 
with the spin excitations and/or disorder. Among all of 
them, the spin-phonon interaction seem to be the best 
candidate. 



The occurrence of a strong Peierls-type Fermi-surface 
instability involving breathing-type displacements of the 
oxygen atoms of the Cu02 planes suggests that spin- 
phonon interactions could contribute to the damping of 
the spin excitations of the insulating compounds. The in- 
plane phonon modes have small energy compared to the 
exchange interaction. Indeed, infrared data for La2Cu04 
gives frequencies of the O stretch modes as 550 and 
690 cm -1 . This corresponds to a bond-stretching force 
constant K of ~ 7.5 and ~ 11 eV/A respectively 0] 
i.e.yj < x 2 > ~ 0.1A The electron-phonon interaction on 
the real material, modifies the exchange interaction pa- 
rameter of the effective spin Hamiltonian mainly through 
the dependence of t p d and the charge-transfer energy on 
the Cu — O distance dcu-o @|- In the spin-wave ap- 
proximation, the lowest order contribution to the damp- 
ing of spin excitations due to the spin-phonon interac- 
tion is proportional to |VJ| 2 [Q. Because along the 
M2CUO4 family, J changes almost linearly with dcu-o, 
|VJ| takes the constant value ~ 4350 cm -1 / A, so the 
maximum contribution to the exchange constant due to 
the spin-phonon coupling is A J ~ ±435 cm -1 . 

The spin-phonon interaction produces, through quan- 
tum fluctuations, some kind of dynamically induced dis- 
order which can be seen by short-wavelength spin exci- 
tations. While the spin dynamic is driven by Hlu s = J, 
the phonon movement is by fru> p h ~ (20 — 60)meV < J, 
then the spin system looks the oxygen atoms as frozen 
in different positions during each spin exchange process. 
Since the phonon energy is much smaller than the ex- 
change parameter, as a hrst approximation, the phonon 
energy can be taken as zero, then the spin-phonon in- 
teraction shows up through fluctuations of the phonon 
field. In this limit phonons can be modeled by static dis- 
order in other words, the fluctuations of the phonon 
field can be modeled by a random distribution P( J, D) 
of the exchange parameter of width D centered around 
J. For the one-band Hubbard model this correspond to 
a distribution for hoppings, centered around t (t') with 
dispersion Dt (Dt 1 ). In what follows we assume a Gaus- 
sian distribution for them. The value of D is of order 
AV < x 2 >/t ~ (0.07 — 0.7)eV for an electron-phonon 
coupling A ~ (0.3 — 3.0) pi) . Since the real value of 
D depends on the details of the microscopic model used 
for the planes as well as the material itself, we fixed 
D = 0.13 which gives the better overall description of 
the mentioned insulators. 

Fig. 2 shows a comparison between experimental and 
calculated Ai g and B\ g spectra for seven different com- 
pounds. Table II presents the calculated and experimen- 
tal values for the first two cumulants. All the spectra of 
Fig. 2 results from a quenched averaged over 150 samples. 
The calculated B\ g response shows many of the exper- 
imentally observed features. Indeed, the characteristic 
2 — magnon peak is located at uj ~ 3 J, the spectra ex- 
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tend up to 8 J and significant spectral weight is found at 
oj ~ (3.5 — 4)J. The A\ g scattering appears in the 
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FIG. 2. Big and Ai g Raman spectra for the 20 sites clus- 
ter. The value of the parameters are given in Table I and 
D = 0.13. Scattered symbols correspond to experimental 
data taken form Ref. ]l2]-]l5t while solid lines are the result of 
theoretical calculations using the effective spin model taking 
into account fluctuations of the phonon field. 

same frequency range of the Bi g response 2 J < uj < 4J 
and the center of gravity of this line is slightly shifted 
to the right respect to the 2 — magnon peak. As it was 
stated in Ref. |Q, the B\ g line can be well fitted with 
two gaussian curves. One centered at w ~ 3J, and the 
other at u> ~ 4 J with a relative intensity of 25% respect 
to the first one. Hence, from the explanation given above 
and the results presented in Ref. j3l) , it is clear that the 
Heisenberg model, even with the inclusion of spin-phonon 



interaction, cannot reproduce simultaneously the width 
and the asymmetry of the B\ g peak. It is necessary to 
take unrealistic values of disorder {D — 0.5J) |3l|] in 
order to transfer spectral weight from the 2 — magnon 
peak to higher energies. From Figs.l and 2, it becomes 
clear that the inclusion of additional terms (J' and K) 
are necessary to reproduce the whole shape of the B lg 
line. These terms generate the whole structure of the A lg 
line because Oa i does not commute with the Heisenberg 
Hamiltonian. This structure has two main contributions 
of approximately same weight. The first is the same as 
the 2 — magnon Bi g peak, now allowed in this repre- 
sentation due to the broken symmetry introduced by the 
fluctuations of the phonon field. Therefore the intensity 
of this peak increases with D. The second contribution 
corresponds to 4 — magnon states (see Fig. 1(b)) mainly 
introduced by the J' and K terms of the Hamiltonian. 
Contrary to the B\ g contribution, this one does not have 
a regular behavior with D. If we apart from D = 0.13 
(the dispersion used to describe all B lg lines) the two 
main A\ g peaks will not have the same intensity and the 
agreement with experiment will pair off. This point is 
clearly seen by comparing Fig. 2 with Fig.l of Ref. ]3l| ] 
where the A\ g line, using a Heisenberg model plus disor- 
der, is given. 

IV. SUMMARY 

Summarizing, we have studied the role of spin-phonon 
and multispin interactions on the line shapes of light- 
scattering experiments within the Loudon-Fleury theory 
for the insulating cuprates. For the insulating parent ma- 
terials of high-temperature superconductors, by contrast 
to conventional antiferromagnets, the energy of short- 
wavelength spin-excitations is greater than the character- 
istic energy of optical phonons so a non-negligible contri- 
bution of the spin-phonon interaction is expected for the 
damping of spin-excitations. Good qualitative agreement 
for the B\ g an the A\ g was found. 

As it was noted before, to describe the dependence of 
the spectrum with the incident laser, a resonant theory 
should be used. But, as it is shown in a recent experiment 
on PrBa2Cu2.7Alo.3O7 [[l6| (see figure 2), resonant effects 
become weak as the laser frequency is decreased. For low 
enough frequencies, the system is out of resonance. The 
main characteristics of the B\ g line measured in this con- 
ditions, the width of the main peak and the shoulder at 
higher energies, remain unaltered. By including spin- 
phonon interactions in the framework of the single-band 
Hubbard model we can reproduce quantitatively these 
features. The A\ g spectra are also well described within 
this theory. Up to our knowledge, this is the first time in 
which a theory with realistic parameters can reproduce 
the Bi g and A\ g measured spectra. Therefore we con- 
clude that the Loudon-Fleury theory is a good starting 
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point to describe the Raman spectra of the insulating 
cuprates. 

In a recent work, a numerical calculation of the Heisen- 
berg model in clusters as big as 6x6 was performed and 
two 2-magnon peaks were found. If this is confirmed for 
bigger clusters, note that Monte Carlo-Max Entropy cal- 
culations on bigger clusters do not show such structures, 
this could lower the phonon-magnon coupling needed to 
adjust the experiments. Unfortunately this is far from 
nowadays numerical possibilities. 

Finite-cluster calculations for 8 insulating cuprates 
support the view of two contributing mechanism to the 
Big and A\ g Raman lines. Although they contribute to 
different parts of the spectra, together, they account for 
the enhanced linewidth and the asymmetry of the B ig 
mode, as well as for a non-negligible response for the A\ g 
Raman line observed in these materials. 
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Table I: Values of the single-band parameters for the different insulators. 





YBa 2 Cu 3 6 . 2 


Pr 2 Cu0 4 


Nd 2 Cu0 4 


Sm 2 Cu04 


La 2 Cu0 4 


Gd 2 Cu0 4 


BSCYCO 


PRBACALO 


u/t 


9.5 


9.5 


9.5 


9.5 


10.5 


9.5 


10.5 


9.5 


t'/t 


0.28 


0.30 


0.28 


0.30 


0.30 


0.30 


0.30 


0.30 


t (eV) 


0.40 


0.42 


0.42 


0.435 


0.475 


0.435 


0.45 


0.37 


J(eV) 


0.138 


0.145 


0.145 


0.150 


0.154 


0.150 


0.146 


0.140 



Table II: Values of the cumulants Mi and M 2 for Raman lines B\ g 
values taken for Ref. [Ol-O. 



and A\ g . First lines correspond to experimental 





YBaCuO 


Pr 2 Cu0 4 


Nd 2 Cu0 4 


Sm 2 Cu0 4 


La 2 Cu0 4 


GdCuO 


BSCYCO 


PRBACAL 




Mi 


1.11 


1.01 


1.03 


1.00 


0.96 


0.97 


0.95 


0.93 






1.02 


0.96 


1.03 


1.00 


0.97 


0.99 


0.93 


0.96 


Bi 


M 2 


0.40 


0.33 


0.37 


0.37 


0.30 


0.32 


0.27 


0.28 






0.42 


0.30 


0.39 


0.40 


0.38 


0.42 


0.27 


0.30 




M 2 /Mi 


0.36 


0.33 


0.36 


0.37 


0.31 


0.32 


0.29 


0.30 






0.41 


0.31 


0.39 


0.40 


0.39 


0.42 


0.29 


0.31 




Mi 


1.30 


1.13 


1.15 


1.04 


1.05 


1.23 


1.20 








1.10 


1.07 


1.13 


1.09 


1.03 


1.10 


1.03 






M 2 


0.42 


0.39 


0.45 


0.37 


0.38 


0.42 


0.46 








0.35 


0.33 


0.37 


0.36 


0.33 


0.37 


0.35 






M 2 /Mi 


0.32 


0.35 


0.39 


0.36 


0.36 


0.34 


0.39 








0.32 


0.30 


0.33 


0.33 


0.32 


0.33 


0.34 
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